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The beta function of the multichannel Kondo model is calculated exactly in the limit of large spin 
00 ' N and channel number M = "fN , with constant 7. There are no corrections in any finite order of 

ON . 1/N. One zero is found at a finite coupling strength, showing directly the Non-Fermi liquid behavior 

ON 1 of the model. This renormalization group flow allows to introduce a variational principle for the 

entropy, to obtain the low temperature thermodynamics. Such in particular the low temperature 
thermodynamics of the non-crossing approximation to the Kondo model becomes accessible. 
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Although the concept of a renormalization group flow in quantum field theory has for a long time been established, 
the existence of a nontrivial fixed point can be shown usually jl| only if it occurs near the trivial one of the nonin- 
teracting theory. It is therefore perhaps worth mentioning that the well studied multichannel Kondo model is such 
a rare case, where the renormalization group flow can be obtained to all orders in the coupling constant, beyond the 
perturbative regime. 

This model describes M identical structureless conduction bands of width D, of fcrmions with an N fold degenerate 
spin degree of freedom, which scatter via an spin-exchange interaction of strength g/N at one localized impurity. It 
is known from a physical argument Q and a 1 /M expansion ^,|| that the renormalization group flow of this model 
should have at least one low temperature stable fixed point at a finite scattering strength, in contrast to the one 
1 1 channel case Q| where it is at infinite coupling. The critical behavior at low temperatures has been calculated for 
some quantities |j|-[3j, showing that the system is a Non-Fermi liquid for M, N > 2, and the critical exponents being 
a function of 7 = M/N only. 

In this work, the result of Refs. is generalized and the (3 function is calculated in the limit of large N with 
7 = M/N kept constant. There will be exactly one intermediate coupling fixed point g c ~ 7 _1 , and (3 1 (g c ) is evaluated 
to be 1/(1 + 7). It is shown explicitly how the renormalization group flow depends on the cutoff scheme between its 
^ fixed points. There will be no corrections to the result in any finite order of 'UN . A variational principle for the low 
. temperature entropy is introduced, generalizing the one for T = in Ref. |7|], and leading to the low temperature 
■^j- thermodynamics. 

■ The multichannel Kondo Hamiltonian is given by 

H= S Tec + c~— ^Tf+cc+f (I) 

On ' 

Here c^ mn creates a conduction electron with spin n = 1, . . . , N and scattering channel m = 1, . . . , M, momentum 
p and energy e p . For the density of states, a sharp cutoff at ±D is used. The fermionic operator /+ describes the 
magnetic impurity configuration \n) = /+|vac). The impurity side is supposed to be singly occupied: fnfn = 1- 
The Boltzmann constant and the density of states at the Fermi surface e p = are set to unity. For g > 0, the 
q ' interaction is antiferromagnetic. 

The diagrammatic technique for this model is standard [0. Each interaction vertex in the scattering channel m is 
, represented by a bare pseudo propagator Sm (z) = g without dynamics, which is for short called a propagator. The 



magnetic configurations have the bare propagators R^\z) = 1/z. 



The full propagators determine completely the physics of H , because every observable can be obtained from their 
Luttinger-Ward functional by coupling to an auxiliary field dJ§]. To calculate the f3 function of H, it is hence enough 
to show that in the low energy limit z, T D, the propagators scale as functions of z, T with the Kondo temperature 
Tk, which will be calculated, depending on the cutoff scheme, leading to a renormalization group flow, with a /3 
function defined as 

chnoTK = -p(g)d g T K . (2) 

It was observed for the related multichannel Anderson model Q that the zero-temperature critical exponents of certain 
observables allow for a controlled expansion in the limit of large M, N with fixed 7. That limit will be introduced 
here as well. It amounts to summing all diagrams with non-crossing conduction electron lines (NCA pj8|), because 
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the 1 /iV-order of a diagram for a propagator is not changed, if a self-energy of second order is eliminated, and just for 
these diagrams this procedure stops if the zero order contribution is reached. Corrections are of relative order 1/N 2 
or smaller B; these are dealt with later. 

The Dyson equations in this limit are obtained by keeping the self-energies of second order in the Luttinger-Ward 
functional but inserting the dressed propagators M 

,D 

S- 1 (z)=g- 1 + f(e)R(z + e)de 

J-D 

r D 

R~ 1 {z) = z + E + 1 f{e)S(z + e)de (3) 

J-D 

where / is the Fermi function, the indices m, n have been dropped, and the ground state energy Eq has been subtracted. 
Scaling for R and S will be established first for frequencies T <C —z <C D and then generalized to all T, \z\ <C D. To 
this end, Eqs. (||) are differentiated with respect to z, neglecting terms oc l/D and 0(exp[— D/T)), 

/oo 
de 5 T (z - e)R(e) (4) 
-OO 

d z R-\z) = 1 - 1 S{z - D) +7 / de S T (z - e)S(e) 



where 5t{z) = —d z f(z). For T <C — z, Eqs. (Q) can be converted into the temperature independent differential 
equations by neglecting terms of the order 0(T/z), 

d z S- 1 (z) = R(z) (5) 
d z R-\z) = 1 - jS(z — D) + jS(z). 

Because Eqs. are independent of g, its value will be fixed implicitly by imposing the boundary condition 

S(-D) d =g. (6) 

g will play the role of a dressed coupling, and coincides in the limit of small coupling with g because of Eqs. (^|). This 
is the analogue of choosing the renormalization point in the usual, perturbative renormalization scheme. The second 
boundary condition is R^ 1 {z = T = 0) = 0, because R^ 1 vanishes at the N — > oo ground state energy Eq which was 
fixed in Eqs. (||) at z — 0. The two boundary conditions determine the solutions of Eqs. (§),(|^). 
Now choose A with r«4«fl. The function 

I(z) = ( 7 - 1 -g)(S- 1 (z)~g~ 1 ) (7) 
+ lniS-^z)^ - 7- 1 \n(R- 1 {z)/R{-D)) 

is because of d z I(z) = R(z)(S(z — D) — g) and Eq. (||) an integral for Eqs. (^) for T <C — z < A if terms of the order 
l/D are neglected. Therefore in this frequency region the boundary condition Eq. (||) can be replaced by the condition 
/ = I(—A). The propagators R, S determine the single particle excitation spectrum of H in the limit N — > 00 [0. 
Hence the solutions of Eqs. (|^) are real and have no singularity below the ground state energy corresponding to z = 0. 
Therefore and because the solutions of Eqs. (|^) depend for z < smoothly on the original boundary condition Eq. (|^), 

the quantity D = f — D) cxp[— A)] is a smooth, real function of g. The rescaling z — > z/D, R^ 1 — > R~ 1 /D 

in Eqs. (||) shows that D/D is a function of g only. If the Kondo temperature is now defined as 

T K = Dg^j- 1 - 5r (7+1) exp[ 7 - 1/g], (8) 

then after the rescaling u = z/Tk and 

r(u) = (7- 1 - 9)T K R(z) (9) 
^) = (7- 1 -ff)- 1 ^) 

the boundary condition / = 7(— A) gives with Eq. (Q) the scale invariant integral 

[-r- 1 (M)] 1/7 = s- 1 (u)exp[s- 1 (u)]. (10) 
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Applying now Eq. (^|) to Eqs. (||), it is evident that it scales with Tk for \z\ <C D, because its integral Eq. ( |To| ) does. 
However, Eqs. (JsJ) is the high frequency limit — z ^> T of Eqs. (||). If therefore additionally T — » T/Tk, its solution 
also scales with Tr- if T, |z| <C £>; hence the solution of Eqs. (|). This shows that there is a renormalization group 
flow, with a low temperature fixed point according to Eq. (B) at g c = 7" 1 . The value <?(g c ) depends via Eq. (^) on 
the particular cutoff-scheme. In terms of the /3 function, defined by Eq. (Q), the main result reads as 

[o ( ~\ = g 2 (7 _1 ~g) nn 

PW 7- 1 +.9 2 + a(ff)5 2 ( 7 - 1 -.9) 1 j 

where a(</) = dghiD/D is a smooth, non universal function of g. This result is still valid when the coupling g c 
becomes large. The then arising nonuniversal terms of the (3 function can be entirely expressed in terms of D and g. 
Two limiting cases are considered in the sequel. For small 7 , assuming g < g ci the expansion in 7 _1 is possible, 

(3(g) = —g 2 + jg 3 + ^yg 4 — j 2 g 5 + 0(7~ 4 ) + nonuniversal terms 

and coincides || with the large M expansion for N — 2. In the other limit g w g c , it is customary to use the running 
coupling constant g(T), defined by TdrgiT) = (3(g(T)) with the initial condition g(D) — g. Its low temperature form 
is determined by the universal number (3'(g c ) = 1/(7 + 1), as 

| 5c -< 7 (T)|ocT 1 /(i+'r). (12) 

Next, the corrections at finite 1/N are considered. The Luttinger-Ward functional is expanded in 1/N, yielding 
diagrams for products of terms R(z)'E v (R(z) , S(z)) with v > 2, where T,^ is the ^th order self-energy of the 
magnetic configuration propagator. In each such diagram, the number of propagators R equals those of S and the 
number of integrations J_ D f(e)de over the energies of the conduction electrons. After the rescaling of T, z, e, R, and 
S with Tk, in the low energy limit z, T <c D as in Eq. (0), all factors Tk and 7 -1 — g are canceling each other. Hence 
the Luttinger-Ward functional still scales with Tk at low energies, so that there are no corrections to the (3 function 
at any finite order in 1/N. 

In order to extract the low temperature thermodynamics of the model, consider its entropy S, first at small coupling 
j<1 and temperatures Tk < T <C D. Then perturbation theory in g can be applied, for each order in T/D JlCp , 

S(T, g) = S (g, T) + S^g, T)T/D + 0(T 2 /D 2 ). (13) 

With Eqs. ©,@ this can be renormalized in terms of the running coupling constant g(T), to give near g c 

S = S (g(T)) + Axl.9, - g{T)\ 1+ ^ + ..., (14) 

where Ai is a constant. Now, if T is lowered the entropy will decrease. But because of the renormalization group flow, 
S will decrease also if g(T) — ■> g c . Therefore S has a minimum at g = g c . This is the generalization of the variational 
principle for the T = entropy . 

There are two cases to be distinguished: In the case 7 > 1 it can be assumed that perturbation theory is valid for 
all g < g c . Then So is dominating, because applying the variational principle to S s» So, it can be expanded around 
its minimum g c to second order, so that with Eq. (|12| ) follows (A2 being a constant) 

S* = 5'o(T = 0) + A 2 T 2 /( 1 +^ + ... . 

This justifies the assumption that So dominates the second term cx T in Eq. (Q). For 7 < 1, the situation is opposite, 
giving for 7 ^ 1 the specific heat 

C{T) oc y min {2/(i+7),i} + , . . ( 15 ) 

In Ref. [0, this law was obtained by a large N expansion of Hamiltonian ([!]) but with the different constraint 
J2 n fnfn — 9oA with constant go, instead of the constraint J2 n fnfn = 1 used here. Although those results cannot 
be expanded in go = the low temperature thermodynamics of the two models coincide. The reason is that the 
variational principle introduced above, works quite generally if the model has a renormalization group flow with a 
finite, low temperature fixed point g c . The low temperature thermodynamics depends then only on I3'(g c ), which 
coincides with the value obtained here. 

To conclude, in this work the renormalization group flow and the [3 function for the multichannel Kondo model have 
been derived, by a large N expansion to all orders in 1/N. The dependence of the (3 function on the cutoff, could be 
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expressed in terms of two smooth functions of the coupling constant. In contrast to the usual perturbative approaches 
to scaling, this method remains valid even if the running coupling constant becomes large. The zero temperature 
variational principle of Ref. |Q was generalized to the leading low temperature expansion of the entropy, which was 
shown to have a minimum at the low temperature fixed point of the coupling constant flow. This in turn yielded the 
low temperature thermodynamics of the model, to all orders in 1/N. This means in particular that Eq. ( |l5| ) describes 
the low temperature thermodynamics of the non-crossing approximation to the SU(N) Kondo model. 
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expresses his gratitude for the financial support by the Japan Society for the Promotion of Science. 
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